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PREFACE 


This  study  was  authorized  as  a  part  of  the  Civil  Works  Research  and 
Development  Program  by  the  Office,  Chief  of  Engineers  (OCE) ,  US  Army.  The 
work  was  performed  under  the  work  unit  C31551,  Numerical  Modeling  of  Shoreline 
Response  to  Coastal  Structures,  which  is  part  of  the  Shore  Protection  and 
Restoration  Program.  Mr.  J.  H.  Lockhart,  Jr.,  was  the  OCE  Technical  Monitor. 

This  study  was  conducted  from  1  October  1982  through  30  September  1983  by 
the  Coastal  Engineering  Research  Center  (CERC)  in  conjunction  with  related 
scientific  studies  by  Professors  Robert  G.  Dean,  University  of  Florida,  and 
Robert  A.  Dalrymple,  University  of  Delaware.  This  report  represents  a  tech¬ 
nical  paper  prepared  to  present  the  overall  results  of  these  efforts.  The 
CERC  portion  of  the  study  was  under  the  general  supervision  of  Dr.  Robert  W. 
Whalin,  former  Technical  Director  and  present  Chief  of  CERC;  Dr.  Lewis  E.  Link, 
Jr.,  Assistant  Chief  of  CERC;  Mr.  Rudolph  P.  Savage,  former  Chief  of  the 
Research  Division;  Dr.  James  R.  Houston,  present  Chief  of  the  Research  Division 
and  Manager,  Shore  Protection  and  Restoration  Program;  Dr.  Craig  T.  Everts  and 
Mr.  H.  Lee  Butler,  former  and  present  Chiefs  of  the  Coastal  Processes  Branch; 
and  Mr.  William  R.  Dally,  Hydraulic  Engineer,  and  Dr.  S.  Rao  Vemulakonda, 
Research  Hydraulic  Engineer,  former  and  present  Principal  Investigators  of  the 
Numerical  Modeling  of  Shoreline  Response  to  Coastal  Structures  work  unit. 
Professor  Kiyoshi  Horikawa  of  the  University  of  Tokyo  provided  the  laboratory 
data  used  in  this  work.  Mr.  Philip  Vitale  generated  the  three-dimensional  and 
contour  plots. 

Commander  and  Director  of  WES  during  the  conduct  of  this  investigation 
and  the  preparation  of  this  report  was  COL  Tilford  C.  Creel,  CE.  Technical 
Director  was  Mr.  F.  R.  Brown. 
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MODELING  WAVE  TRANSFORMATION  IN  THE  SURF  ZONE 


PART  I :  INTRODUCTION 

Statement  of  t!  e  Problem 

1.  A  major  problem  encountered  in  modeling  nearshore  wave-induced  phe¬ 
nomena  is  the  description  of  wave  parameters  subsequent  to  the  initiation  of 
wave  breaking.  Specifically,  wave  height  and  its  spatial  gradients  generate 
or  have  direct  impact  on  sediment  mobilization  and  suspension,  littoral  cur¬ 
rents  (in  both  the  alongshore  and  on/offshore  directions),  wave-induced  setup, 
and  forces  on  coastal  structures.  While  the  0.78  criterion  (ratio  of  breaker 
height  to  water  depth  =  0.78)  appears  to  provide  a  reasonable  prediction  of 
incipient  breaking  on  mildly  sloping  beaches,  data  from  several  investigators 
(Horikawa  and  Kuo  1966,  Nakamura,  Shiraishi,  and  Sasaki  1966,  Street  and 
Camfield  1966,  Divoky,  Le  Mehaute/,  and  Lin  1970)  show  that  this  criterion  does 
not  hold  farther  into  the  surf  zone  and  that  such  a  similarity  model  is  more 
inappropriate  cm  mild  slopes  -  just  where  many  coastal  scientists  assume  it  is 
most  valid.  Another  shortcoming  of  this  and  many  other  representations  de¬ 
veloped  to  date  is  that  they  are  not  applicable  on  non-monotonic  beach  pro¬ 
files  such  as  those  containing  bar/trough  formations. 

Purpose  of  the  Study 


2.  The  purpose  of  this  study  is  to  develop  a  general  numerical  model  for 
wave  decay  and  transformation  subsequent  to  breaking  which  includes  the  ef¬ 
fects  of  bottom  friction,  wave-induced  setup,  and  beach  profiles  of  arbitrary 
shape.  The  model  will  be  calibrated  and  verified  with  laboratory  and  proto¬ 
type  scale  data. 


PART  II:  MODEL  DEVELOPMENT 


Literature  Reviev 


3.  During  the  past  two  decades  a  number  of  laboratory,  field,  and  analy¬ 
tical  studies  have  been  conducted  to  develop  a  realistic  model  of  wave  height 
transformation  across  the  surf  zone.  Such  a  model  is  essential  to  an  adequate 
understanding  of  nearshore  hydrodynamics  and  longshore  sediment  transport. 

The  steady-state  equation  governing  energy  balance  for  waves  advancing 
directly  toward  shore  is* 

3(EC  ) 

-  —  6  (x )  (1) 


in  which 

E  =  wave  energy  per  unit  surface  area 


g 

x 


group  velocity 

horizontal  coordinate  normal  to  the  beach  and  increasing  toward  the 
shore 


6  =  energy  dissipation  rate  per  unit  surface  area  due  to  boundary  shear, 
turbulence  due  to  breaking,  etc. 

The  central  problem  in  previous  studies  has  been  the  development  of  a  rational 
and  universally  valid  formulation  for  6  .  The  most  physically  appealing  ap¬ 
proach,  first  advanced  by  Le  Mahanfe'  (1962),  has  been  the  approximation  of  a 
breaking  wave  as  a  propagating  bore  (or  hydraulic  jump)  in  which  case  6  is 
given  by 


6  -  p 
4  ,  2 


(2) 


where 

o  =  mass  density  of  water 
g  =  acceleration  due  to  gravity 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  B) . 


B  =  parameter  representing  the  fraction  of  the  wave  height  that  is  due  to 
breaking 

H  =  wave  height 

h  =  water  depth 

Q  =  transport  of  water  across  the  bore 
Since  its  first  introduction  by  Le  Mehaute,>  Equation  2  has  been  used  in 
slightly  different  forms  by  a  number  of  investigators  to  represent  periodic 
water  waves  in  the  laboratory  (Divoky,  Le  Mehaute\  and  Lin  (1970),  Hwang  and 
Divoky  (1970),  Svendsen,  Madsen  and  Hansen  (1978)),  aperiodic  water  waves  in 
the  laboratory  (Battjes  and  Janssen  (1978)),  and  aperiodic  waves  in  nature 
(Thornton  and  Guza  (1983)). 

4.  Other  approaches  for  wave  energy  dissipation  have  included  that  of 
Horikawa  and  Kuo  (1966)  in  which  the  internal  energy  dissipation  is  repre¬ 
sented  in  terms  of  turbulent  velocity  fluctuations  which  are  assumed  to  decay 
exponentially  with  distance  from  the  wave  break  point.  In  a  different  ap¬ 
proach,  Mizuguchi  (1980)  applies  the  analytical  solution  of  Lamb  (1932)  for 
internal  energy  dissipation  due  to  viscosity,  replacing  the  molecular  kine¬ 
matic  viscosity  with  eddy  viscosity  which  must  be  estimated  based  on  the  wave 
and  beach  profile  characteristics. 

5.  Understandably,  the  extension  of  approaches  developed  for  periodic 
waves  to  aperiodic  waves  introduces  complexities,  primarily  with  respect  to 
representation  of  the  probability  distribution  function.  In  a  laboratory 
study  of  breaking  random  waves  over  plane  and  barred  beaches,  Battjes  and 
Janssen  (1978)  employed  Equation  2  interpreting  H  as  H  and  introduced 
two  assignable  constants.  A  truncated  Rayleigh  wave  height  distribution  was 
assumed  with  a  finite  probability  of  the  maximum  truncating  wave  height  oc¬ 
curring,  resulting  in  a  delta  function  at  this  limit.  Battjes  and  Janssen 
included  the  effect  of  wave  setup  and  demonstrated  good  agreement  between  the 
predicted  and  the  measured  wave  height  distribution.  Following  Collins 
(1970),  Kuo  and  Kuo  (1974)  modified  the  truncated  wave  height  probability  dis¬ 
tribution  by  omitting  the  della  function  at  the  upper  wave  height.  Goda 
(1975)  adopted  a  modified  probability  distribution  in  which  the  distribution 
decreases  linearly  to  zero  over  a  specified  upper  range  of  the  wave  heights. 
Based  on  field  measurements,  Thornton  and  Guza  (1983)  present  convincing  data 
that  the  unmodified  Rayleigh  distribution  is  applicable  seaward  of  and  across 


the  entire  surf  zone.  At  each  location,  they  define  a  breaking  subpopulation 
of  the  waves. 

6.  In  the  course  of  the  investigations  described  above,  a  substantial 
data  base  of  wave  height  transformation  across  laboratory  and  field  beaches  of 
various  profiles  has  been  developed.  Nakamura,  Shiraishi  and  Sasaki  (1966) 
presented  data  based  over  1500  runs  of  wave  transformation  across  a  plane 
beach  in  the  laboratory.  This  assemblage  of  wave  data  is  useful  in  providing 
guidance  and  a  basis  for  evaluation  in  an  attempt  to  develop  a  universally 
valid  model. 

7.  None  of  the  models  developed  and  evaluated  to  date  provide  a  demon¬ 
strated,  completely  general  capability  for  representing  wave  transformation 
across  the  surf  zone.  The  waves  in  this  zone  are  not  linear  and  some  of  the 
difficulty  may  be  in  the  application  of  linear  theory,  although  other  theories 
(cnoidal  and  solitary)  have  been  used  by  some  investigators. 

8.  One  of  the  features  not  represented  in  most  models  is  that  of  the 
wave  height  stabilizing  at  some  value  in  a  uniform  depth  following  the  initia¬ 
tion  of  breaking.  The  laboratory  data  of  Horikawa  and  Kuo  (1966)  and  general 
observations  and  intuition  support  such  a  phenomenon,  yet  none  of  the  energy 
dissipation  models  based  on  the  moving  hydraulic  jump  predicts  this  effect. 
Although  the  model  by  Mizuguchi  (1980)  includes  this  stabilization,  a  simpler 
model  could  minimize  the  difficulties  in  estimating  the  representative  eddy 
viscosity  and  in  rationally  applying  the  model  to  a  beach  of  nonuniform  slope. 

9.  This  report  centers  on  the  development  and  evaluation  of  a 
model  which  includes  the  wave  height  stabilization  and  which  appears  to  be 
broadly  applicable  based  on  comparison  with  periodic  water  wave  data  from 
laboratory  studies.  In  comparison  with  the  laboratory  data,  the  available 
field  data  are  much  more  limited.  The  evaluation  of  a  model  by  laboratory 
data  would  further  an  understanding  of  the  problem. 


Theoretical  Development 


10.  The  primary  purpose  of  this  work  is  to  accurately  describe  the  decay 
in  wave  height  due  to  breaking  as  a  wave  crosses  the  surf  zone.  Ultimately, 
it  is  desirable  to  address  as  much  of  the  wave  transformation  process  as  pos¬ 
sible,  including  shoaling,  breaking,  and  reformation.  Due  to  the  complexity 
of  the  problem,  a  common  approach  is  to  draw  an  analogy  between  the  turbulent 


V.VA.V/.  A  *>  * 


front  face  of  a  breaking  wave  and  a  moving  hydraulic  jump  (Le  Mehaute*  (1962), 
Peregrine  and  Svendsen  (1978)).  However,  in  the  present  model  a  more  macro¬ 
scopic  treatment  is  applied  where  the  energy  dissipation  across  the  entire 
surf  zone  is  compared  to  that  of  a  hydraulic  jump.  Rather  than  attempting  to 
quantify  the  local  energy  dissipation  by  addressing  the  turbulent  velocity 
fluctuations,  a  common  method  is  to  describe  the  local  gradient  in  depth- 
integrated  energy  flux  as  stated  in  Equation  1.  In  the  typical  hydraulic  jump 
in  a  rectangular  channel  (Figure  la),  energy  is  dissipated  over  some  distance 
to  transform  a  high  energy  condition  (Froude  Number  fF  >  1)  to  a  lower  energy 
condition  (  fF  <  1)  by  way  of  turbulence.  The  specific  energy  (energy  per  unit 
weight)  of  the  flow  anywhere  along  the  channel  as  introduced  by  Bakhmeteff  in 
1912  is  given  by 


2 

E'  =  h  +  (3) 

where  it  is  assumed  that  the  horizontal  fluid  velocity  u  is  uniform  over 
the  depth  of  flow  h  .  The  energy  per  unit  volume  is  then 

2 

E  =  yh  +  (4) 

where  y  is  the  unit  weight.  The  energy  flux  per  unit  width  of  the  channel 
past  any  section  is  given  by 


"Flux 


Euh 


yuh2  + 


(5) 


Defining  the  flow  rate  per  unit  width  of  channel,  which  is  spatially  and 
temporally  constant,  as  q  =  uh  ,  Equation  5  becomes 


E  - 
Flux 


yqh  + 


(6) 


The  local  rate  of  energy  dissipation  per  unit  plan  area  is  simply  the  spatial 
gradient  of  the  energy  flux 
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dE 


Flux 

dx 


h  pgM  2  dh 
2  *  2h2/hdx 


(7) 


which  can  also  be  written  as 


dE 


Flux 

dx 


2  dh 
h  dx 


(8) 


The  term  in  parentheses  in  Equation  8  is  the  local  energy  flux  as  previously 
defined  by  Equation  6.  The  second  term  in  the  brackets  is  equivalent  to  the 
energy  flux  that  would  be  present  if  the  flow  were  critical  (  F  =  1)  at  that 
location.  Moreover,  if  the  flow  were  critical,  this  would  be  the  minimum 
energy  flux  that  would  occur  for  the  given  flow  rate.  Equation  8  is  then 


H 

The  local  slope  of  the  water  surface  dh/dx  cannot  be  determined  without  a 
knowledge  of  the  details  of  the  vertical  flow  regime  within  the  jump.  For¬ 
tunately,  it  is  the  form  of  Equation  9  that  will  be  of  interest  when  the 
breaking  wave  case  is  addressed. 

11.  Next,  consider  a  beach  profile  that  rises  from  deco  water  in  a 
gently  sloping  manner  and  at  some  point  in  shallow  water  becomes  horizontal 
(see  Figure  lb).  Consider  further  a  wave  propagating  onto  this  profile  with 
characteristics  such  that  breaking  starts  at  the  point  where  the  bottom  be¬ 
comes  horizontal.  The  wave  will  not  instantaneously  stop  breaking  because  the 
bottom  has  become  horizontal  (as  dictated  by  the  0.78  criterion)  but  would 
continue  to  break  until  some  stable  wave  height  is  attained.  Breaking  would 
be  most  intense  just  shoreward  of  line  AA  and  would  decrease  until  an  approx¬ 
imate  stable  wave  height  is  reached  at  line  BB.  This  idealized  surf  zone  can 
be  thought  of  as  an  elongated  hydraulic  jump,  where  again,  a  high  energy  free 
surface  flow  condition  is  converted  to  a  lower  energy  one  by  the  dissipation 
of  energy  through  turbulence.  In  both  situations,  it  appears  that  the  fluid 
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flow  surpasses  some  criterion  and  then  dissipates  energy  until  the  end  result 
is  a  flow  condition  well  below  the  criterion.  (For  the  hydraulic  jump,  this 
criterion  is  F  =  1,  but  if  such  a  criterion  exists  for  breaking  waves,  it  is 
heretofore  unresolved.)  For  these  reasons,  the  following  expression,  similar 
to  Equation  9,  will  be  used  to  describe  energy  dissipation  across  the  surf 


^  ECg  -  (ECg)s 


ECg  =  depth-integrated  time-averaged  energy  flux  as  given  by  shallow 
water  linear  wave  theory 

K  =  dimensionless  decay  coefficient 

h'  =  still-water  depth 

ECg  =  energy  flux  associated  with  the  "stable"  wave  condition  that  the 
breaking  wave  is  striving  to  attain 

Horikawa  and  Kuo  (1966)  conducted  laboratory  tests  with  a  bottom  configuration 
identical  to  the  one  described.  As  shown  in  Figure  2,  their  data  indicate  a 
stable  wave  criterion  given  by 


H  =  IV 
s 


where 

H  =  stable  wave  height 
s 

r  =  dimensionless  coefficient  whose  value  appears  to  lie  somewhere 
between  0.35  and  0.40 

12.  In  Figure  3,  wave  height  was  plotted  against  still-water  depth  for  a 
uniform  beach  slope  of  1/65,  revealing  that  the  breaking  waves  tended  to  ap¬ 
proach  asymptotically  the  line  H  =  0.5h'  where  H  is  the  local  wave  height. 
Equation  11  appears  to  be  a  reasonable  assumption,  and  Equation  10  can  then  be 
written  as 


d[H2 (h’)1/2] 


=  [H2  (h’ )  1/2  -  r2(h’)5/2] 

n 
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•  Data  from  Horikawa  and  Kuo 
(1966) 

.^Analytical  Solution(K=0.2,r=0. 


in  which  m  is  the  beach  slope,  and 


■(§  - :)  V- 


Note  that  this  solution  fails  when  K/m  =  5/2  .  For  this  special  case,  the 


solution  is 


where 


Also  note  that  if  K  is  set  equal  to  zero  (no  decay).  Equation  14  becomes 


1L  Vh  7 


which  is  Green's  Law.  If  a  =  -1.0  ,  Equation  14  reverts  to  the  common  sim¬ 
ilarity  model  H  a  h'  .  Equations  14  and  16  are  plotted  in  Figures  4a  and  4b 
for  several  values  of  H^/h^  and  K/m  .  Figure  5  compares  Equation  14  with 
the  data  presented  in  Horikawa  and  Kuo  (1966) . 


Additional  Factors 


15.  During  initial  examination  of  the  complete  raw  data  set  collected  by 
Horikawa  and  Kuo,  it  was  noticed  that  in  all  cases  where  measurements  were 
taken  in  the  inner  portion  of  the  surf  zone,  as  the  still-water  depth  ap¬ 
proached  zero  the  wave  height  did  not  (Figure  5).  To  model  this  phenomenon 
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Figure  5.  Comparison  of  analytical  solution  (Equation  14)  with  wave 
decay  data  as  presented  in  Horikawa  and  Kuo  (1966)  for  various  beach 

slopes  (K  =  0.17  ,  T  =  0.5) 

better,  wave- induced  setup  and  setdown  of  the  mean  water  level  must  be  in¬ 
cluded.  The  same  conclusion  was  reached  originally  by  Hwang  and  Divoky 
(1970).  From  Longuet-Higgins  and  Stewart  (1963),  the  slope  of  the  mean  water 
level  n  is  given  by 


s%  V.  vs  v 


STB 


dn  _  -1 

dx  _  _ , 


pg(h'  +  n) 


where  the  radiation  stress  S  in  shallow  water  is 

xx 


3  2 

S  =  77  Pgn 
xx  16 


Equation  19  becomes 


dn  _  -3  1  dH 

dx  “  16  (h*  +  n)  dx 


and  can  be  used  in  conjunction  with  a  slightly  different  form  of  Equation  10 


iifSi  -  [ECg  -  (ECg)sJ 


and  Equation  12 


dt"2<h>1/2>  -  4  [H2(h)1/2  -  r2(h)5/2J 

dx  h 


in  which  h  ,  the  mean  water  depth  given  by 


h  =  h'  +  n  (24 

has  replaced  the  still-water  depth  h'  . 

Bottom  friction 

16.  Although  energy  dissipation  due  to  bottom  friction  is  generally 
negligible  when  compared  to  energy  loss  due  to  breaking,  the  former  will  be 
examined  in  this  paragraph.  The  average  rate  of  energy  dissipation  per  plan 
area  due  to  bottom  friction  is  expressed  by 


•  •  »  •  ■  •  •  *  '  *«  *  <  «  J 

.  •  .*  fcS'-'ji  *  *  -r  V  V  .*  A  / 


WWW 


(25) 


ELB=  T  J  TBUBdt 


where  T  is  the  wave  period,  xD  is  the  bottom  shear  stress  given  by 

D 


TB  =  p  2  VV 


(26) 


and  uD  is  taken  to  be  the  first-order  horizontal  water  particle  velocity  at 

D 

the  bottom  which  (from  linear  wave  theory  in  shallow  water)  is 


u 


B 


(27) 


and  t  is  the  time. 

f  =  drag  coefficient  dependent  on  flow  and  bottom/sediment 
characteristics 

o  =  wave  angular  frequency 

Substituting  Equations  26  and  27  into  Equation  25  and  integrating  yields 


e  =  fifsi 

LB  12tt 


(28) 


Numerical  Solution 


17.  Introducing  setup,  beach  profiles  of  more  realistic  shape,  or  bottom 
friction  to  the  model  render  the  equations  unsolvable  analytically.  A  numeri¬ 
cal  scheme  will  therefore  be  developed  which  describes  the  one-dimensional 
transformation  of  wave  height  over  bottoms  of  arbitrary  shape  due  to  shoaling, 
breaking,  reformation,  and  bottom  friction,  including  the  effects  of  setup  in 
mean  water  level. 

18.  Rewriting  Equation  22  in  a  finite  difference  form  using  a  central 
average  for  each  of  the  quantities  on  the  right-hand  side  gives 
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+  (ECg) 


h .  +  n .  +  h 


and  the  superscript  (~)  denotes  conditions  before  bottom  friction  is  applied 
By  applying  shallow-water  wave  theory  and  the  stable  wave  criterion  Hg  =  Th 
Equation  29  may  be  transformed  to 


(1  -  A) 

(1  +  A) 


in  which 


2  2 

F  hi+l 


h5/2 
,2  hi 


Equation  21  is  treated  in  a  similar  manner 


~2  2 

,  H .  .  -  H. 

-  -3  l+l  l  ,  - 

ni+l  16  h  ni 

m 


Before  can  be  calculated,  the  mean  water  level  at  i+1  must  be  known 
The  assumption  is  made  that  the  mean  water  depth  at  the  next  location  can  be 
approximated  by 


i+1 


A  first  estimate  of  can  obtained  using  Equation  31,  after  which  a 

better  estimate  of  h^+^  is  made  using  Equation  33.  The  program  iterates 
until  the  new  value  for  the  mean  water  depth  is  very  close  to  the  old  one.  In 
the  calibration  runs,  no  more  than  three  iterations  were  required  for  the 
difference  in  estimates  to  become  less  than  a  millimeter. 

19.  The  decay  in  wave  energy  due  to  bottom  friction  is  not  included 
directly  in  the  wave  height  computation  but  is  calculated  separately  using 
Equation  28  and  energy  flux  considerations.  This  procedure  yields 


Hi+1  +  Hi 


3,h3/2h^ 

m  i+1 


where  Hi+^  is  the  wave  height  at  the  next  point  of  interest  after  shoaling, 
breaking,  and  bottom  friction  have  been  applied.  Separating  the  bottom  fric¬ 
tion  in  this  manner  is  justified  because  it  plays  a  significant  role  only  in 
the  absence  of  breaking. 

20.  To  apply  the  model  in  a  given  situation,  the  following  information 
is  required:  (a)  the  wave  height  and  still-water  depth  at  a  known  nearshore 
location,  (b)  the  ratio  of  wave  height  to  water  depth  at  incipient  breaking, 
(c)  the  bottom  friction  coefficient,  and  (d)  the  bottom  profile.  The  ratio  of 
breaking  height  to  depth  is  not  easily  predicted  and  was  not  treated  in  an 
extensive  manner  in  this  study.  Mallard  (1978)  offers  a  more  complete  inves¬ 
tigation.  Assuming  the  starting  point  is  in  shallow  water  and  outside  the 
surf  zone,  the  setdown  in  mean  water  level  as  given  by  Longuet-Higgins  and 
Stewart  (1963)  is 


21.  From  these  initial  conditions  and  using  Equations  31,  33,  and  35  as 
described,  the  wave  height  will  increase  (with  some  loss  in  height  due  to  bot¬ 
tom  friction)  as  the  wave  moves  shoreward  until  the  incipient  breaking  cri¬ 
terion  is  reached.  The  wave  then  breaks  until  it  becomes  locally  stable 
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(if  at  all).  On  barred  profiles,  the  combination  of  wave  decay  and  increasing 
water  depth  as  the  wave  passes  over  the  trough  enables  the  wave  to  reach  sta¬ 
bility.  The  "reformed"  wave  then  shoals  again  until  the  breaking  criterion  is 
reached  and  the  process  repeats  until  the  mean  water  depth  reaches  an  arbi¬ 
trarily  chosen  small  value  (0.25  m  is  a  reasonable  choice  at  prototype  scale). 
A  flow  chart  of  the  wave  transformation  model  is  shown  in  Figure  6. 


in 


U 


Calibration  of  the  Model 

22.  The  model  is  calibrated  by  determining  the  best  values  for  the 
stable  wave  factor  (T)  and  the  wave  decay  factor  K  using  a  least-squares 
procedure.  Horikawa  and  Kuo's  (1966)  original  laboratory  data  of  waves 
breaking  on  plane  slopes  were  examined.  Starting  at  incipient  breaking,  the 
location  of  wave  heights  were  measured  by  Horikawa  and  Kuo  at  known  distances 
across  the  surf  zone  under  monochromatic  wave  conditions  for  plane  smooth 
beaches  of  rubber  and  concrete  having  slopes  of  1/20,  1/30,  1/65,  and  1/80. 
The  wave  period  varied  from  1.2  to  2.2  sec  and  the  incipient  breaker  height 
from  7  to  27  cm.  Although  not  specifically  stated  by  Horikawa  and  Kuo,  the 
breakers  must  have  spanned  both  the  plunging  and  spilling  types  because  the 
ratio  of  wave  height  to  water  depth  at  incipient  breaking  ranged  from  0.63  to 
1.67.  The  number  of  waves  and  data  points  for  the  slopes  analyzed  are 
presented  as  follows: 


Slope 

Number  of  Waves 

Number  of  Data  Points 

1/30 

17 

173 

1/65 

12 

96 

1/80 

56 

500 

Data  from  the  1/20  slope  were  not  included  in  the  calibration  because  the  mea¬ 
surements  were  made  too  far  apart  for  the  model  to  remain  numerically  stable. 
The  error  function  to  be  minimized  is  defined  by 


e(I\  K)  = 


£ [H  .(r,  K) 
.1  =  1  PJ 

N  2 
£  % 
i=i 


where 


H  =  measured  wave  height 

H  .  =  wave  height  at  that  location  as  predicted  by  the  numerical  scheme 
for  given  incipient  conditions  and  values  of  T  and  K 

N  =  number  of  data  points  analyzed 


A  .%  N 


A  lengthy  and  involved  attempt  to  obtain  best-fit  values  of  T  and  K  was 
made  using  a  nonlinear,  least-squares  error  analysis.  This  method  involves 
choosing  initial  values  for  the  factors,  fitting  a  parabolic  surface  to  the 
error  surface  at  that  point,  locating  the  minimum  of  the  fitted  surface,  and 
using  these  new  values  of  T  and  K  for  the  next  attempt.  In  this  manner, 
the  procedure  theoretically  converges  on  the  minimum  of  the  error  surface. 
However,  this  method  was  unsuccessful,  apparently  because  the  error  surfaces 
are  too  highly  nonlinear,  i.e.,  a  paraboloid  is  a  poor  approximation  for  the 
error  surface.  By  calculating  the  error  at  regular  intervals  of  K  and  T  , 
a  discretized  error  surface  can  be  generated,  its  low  point  occurring  near  the 
best-  fit  values  for  the  factors.  Figures  7,  8,  and  9  display  three- 
dimensional  and  contour  plots  of  the  error  surfaces  for  the  data  taken  on  the 
1/30,  1/65,  and  1/80  slopes.  Note  the  recurved  shapes  of  all  the  surfaces  and 
the  saddle  point  present  in  two  of  the  contour  plots.  A  one-wave,  sixteen- 
point  artificial  "data"  set  was  generated  using  the  model  on  a  1/80  slope  with 
T  set  equal  to  0.35  and  K  equal  to  0.10  (the  best  fit  values  from  the  real 
1/80  data  set).  Figure  10  displays  the  corresponding  error  surface  and  con¬ 
tour  plots.  Note  that  it  has  the  same  general  shape  and  saddle  point  as  the 
real  data  set.  It  can  therefore  be  concluded  that  the  shape  and  saddle  point 
are  characteristics  of  the  model  and  are  not  due  to  problems  with  the  data 
set.  The  best  fit  values  for  T  and  K  for  the  three  slopes  analyzed  are 
presented  as  follows: 


Slope 

r 

K 

Minimum  Error 

1/30 

0.475 

0.275 

0.1165 

1/65 

0.355 

0.115 

0.1049 

1/80 

0.350 

0.100 

0.1298 

The  best-fit  values  for  the  two  factors  do  vary  with  beach  slope,  especially 
as  the  beach  becomes  steeper.  However,  it  would  be  preferable  to  choose  sin¬ 
gle  values  for  T  and  K  which  give  satisfactory  results  for  all  beach 
slopes,  allowing  the  model  to  be  used  on  beach  profiles  of  more  realistic 
shape.  Fortunately,  the  error  surfaces  for  the  three  slopes  tested  are  rela¬ 
tively  broad  and  flat  in  the  vicinity  of  their  minima,  so  the  factors  can  be 
changed  without  excessively  increasing  the  combined  error.  The  procedure 
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Figure  10.  Three-dimensional  and  contour  plots  of  error  e  between  numerical 
solution  and  model-generated  "data"  as  a  function  of  wave  decay  factor  K  and 
stable  wave  factor  f  for  1/80  beach  slope 


followed  was  to  superimpose  the  contour  plots  of  Figures  7,  8,  and  9  and  find 
the  location  where  the  sum  of  the  three  error  values  is  minimized.  This  point 
occurs  at 

T  =  0.40 

K  =  0.15 

and  it  is  recommended  that  these  values  be  used  in  situations  where  the  bottom 
slope  varies  over  a  wide  range.  If  the  beach  is  nearly  planar,  the  best-fit 
values  presented  for  r  and  K  for  the  three  slopes  may  be  used  accordingly. 

23.  The  breaker  model  was  not  calibrated  to  prototype  conditions  due  to 
a  lack  of  suitable  data.  It  would  require  the  measurement  of  the  breaker 
height  distribution  across  the  surf  zone  under  truly  monochromatic  conditions. 
The  waves  should  be  incident  normal  to  the  bottom  contours,  and  the  beach  pro¬ 
file  must  be  known  and  preferably  monotonic.  A  large  wave  tank  with  a  regular 
wave  generator  would  be  ideal  for  this  type  of  experiment.  Calibration  of  the 
model  in  its  present  form  using  irregular  waves  would  require  following  each 
individual  breaker  as  it  travels  across  the  surf  zone  (provided  wave-wave  in¬ 
teraction  is  not  significant).  Because  the  bubbles  which  dissipate  the  energy 
in  breaking  waves  may  not  be  scaled  down  properly  in  the  laboratory,  it  is  hy¬ 
pothesized  that  K  will  assume  a  lesser  value  under  prototype  conditions.  It 
is  unclear  what  change,  if  any,  would  occur  in  T  . 


PART  III:  RESULTS  FOR  LABORATORY  CONDITIONS 


Wave  Height 

24.  Figures  11-14  display  a  representative  sample  of  model-predicted 
breaker  decay  as  compared  with  the  aforementioned  laboratory  data  for  plane 
beaches  of  slopes  1/20,  1/30,  1/65,  and  1/80.  The  figures  are  dimensional 
plots  of  wave  height  versus  still-water  depth.  In  all  cases,  the  wave  decay 
factor  K  was  set  equal  to  0.15  and  the  stable  wave  factor  r  was  taken  to 
be  0.40.  Bottom  friction  was  considered  negligible.  Each  curve  was  generated 
by  taking  as  input  the  wave  height  and  still-water  depth  at  incipient  breaking 
as  given  in  the  data  and  by  calculating  stepwise  (Ax  =  1/m)*  the  setup  and 
decay  profiles.  It  is  not  practical  to  display  the  data  and  model  results  in 
dimensionless  form  due  to  the  dependency  of  the  shape  of  the  decay  profile  on 
the  incipient  conditions  as  well  as  the  setup.  Moreover,  setup  measurements 
were  not  included  in  the  Horikawa  and  Kuo  experiments. 

25.  Examination  of  these  results  shows  that  the  model  developed  in  this 
study  appears  to  provide  a  good  representation  of  breaking  wave  decay  on  plane 
beaches  of  laboratory  scale.  It  is  important  to  note  that  the  model  is  in 
good  quantitative  agreement  over  the  wide  range  of  slopes  tested,  even  with 
the  two  factors  K  and  T  held  constant  at  values  that  are  not  necessarily 
the  best-fit  values  for  that  particular  slope. 

26.  The  line  H  =  0.78  h'  is  plotted  in  each  of  the  figures  and  appears 
to  be  a  reasonable  description  of  breaker  decay  for  only  the  1/30  slope.  In 
fact,  the  surf  similarity  model  (H  h’)  which  is  so  prevalent  in  the  coastal 
literature  seems  to  be  remotely  valid  only  for  beaches  of  much  greater  slope 
than  those  commonly  found  in  nature. 

27.  For  waves  having  similar  ratios  of  incipient  breaker  height  to  depth 
but  different  incipient  breaker  heights,  the  measured  decay  profiles  inter¬ 
sect.  This  phenomenon  is  apparent  also  in  the  profiles  produced  by  the  model. 
The  greater  total  setup  associated  with  an  initially  larger  wave  (Bowen  1969) 
provides  a  greater  mean  water  depth  and  so  a  larger  stable  wave  height  in  the 
inner  surf  zone.  At  this  point  the  larger  wave  decays  less  rapidly,  causing 
the  decay  profiles  to  intersect. 

*  See  notation  (Appendix  B) . 
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Figure  12.  Comparison  of  breaker  model  and  0.78  criterion  with  the 
laboratory  data  of  Horikawa  and  Kuo  for  1/30  beach  slope 
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Figure  13.  Comparison  of  breaker  model  and  0.78  criterion  with  the 
laboratory  data  of  Horikawa  and  Kuo  for  1/65  beach  slope 
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Figure  14.  Comparison  of  breaker  model  and  0.78  criterion  with  the 
laboratory  data  of  Horikawa  and  Kuo  for  1/80  beach  slope 


28.  Although  the  range  in  wave  period  in  the  data  is  limited  (1.2  - 
2.26  sec),  it  appears  that  wave  period  is  not  a  primary  factor  in  the  decay  of 
wave  height  after  breaking  is  initiated  due  to  the  fact  that  breaking  is  a 
phenomenon  restricted  only  to  a  portion  of  the  wave  crest.  Wave  period  does 
affect  the  wave-height-to-water-depth  ratio  at  incipient  breaking  (along  with 
beach  shape,  slope,  and  deepwater  wave  height)  and  therefore  affects  the  shape 
of  the  decay  profile  through  the  initial  condition.  Wave  period  plays  no 
direct  role  in  the  decay  portion  of  the  present  model  because  the  incipient 
conditions  are  given  as  input  to  the  model. 

29.  Characteristically,  waves  of  the  plunging  type  dissipate  a  large 
portion  of  their  energy  in  a  concentrated  manner  in  the  region  just  shoreward 
of  the  breaker  line,  while  spilling  breakers  dissipate  their  energy  at  a 
slower  rate.  At  incipient  breaking,  the  wave-height-to-depth  ratio  is  usually 
greater  than  1.0  for  plunging  breakers  and  0.78  or  less  for  spilling  breakers. 
Figure  15  presents  predicted  wave  height  decay  curves  for  two  waves  with  the 
same  initial  breaking  depth  but  different  wave  heights.  The  larger  wave  ini¬ 
tially  decays  at  a  much  greater  rate  than  the  smaller,  and  it  appears  the 
model  is  at  least  qualitatively  correct  in  dealing  with  the  breaker  types, 
plunging  and  spilling. 

30.  Figure  16  demonstrates  the  negligible  effect  energy  losses  from  bot¬ 
tom  friction  have  on  the  wave  decay  profile  when  compared  to  breaking.  The 
upper  curve  is  a  test  case  of  the  model  without  bottom  friction.  The  lower 
curve  was  generated  with  the  same  conditions,  except  greatly  exaggerated  bot¬ 
tom  friction  losses  (assuming  no  bed  form  losses)  were  included.  The  drag 
coefficient  f  for  the  bottom  was  set  equal  to  0.10,  about  two  orders  of 
magnitude  greater  than  is  realistic  for  the  smooth  rubber  and  concrete  slopes 
used  by  Horikawa  and  Kuo.  This  was  for  the  purpose  of  distinguishing  between 
the  two  curves. 


31.  The  data  of  Horikawa  and  Kuo  do  not  include  measurements  of  setup, 
so  experimental  data  presented  in  Bowen,  Inman,  and  Simmons  (1968)  and  Bowen 
(1969)  were  examined.  In  these  experiments  wave  heights  and  mean  water  levels 
were  measured  on  a  relatively  steep  plane  beach  of  1/12  slope.  The  results  of 
two  tests  are  presented  in  Figures  17  and  18,  along  with  the  decay  and  setup 


Beach  Slope  =  1/30 


Figure  15.  Comparison  of  model-predicted  wave  decay  for  spilling 

and  plunging  breakers 

curves  predicted  by  the  numerical  model  (K  =  0.25  ,  T  =  0.35).  The  breaker 
decay  compares  well  and  the  maximum  setup  values  are  reasonable  if  swash  zone 
is  neglected;  however,  the  predicted  setup  curves  do  not  follow  the  data. 
Apparently,  Equation  20  is  not  a  good  representation  of  the  onshore  excess 
momentum  flux  for  near-breaking  and  breaking  conditions  as  might  be  expected. 
Higher  order  wave  theories  yield  significantly  less  momentum  flux  for  a  given 
wave  height  than  linear  theory  (Dean  1974),  and  this  difference  is  the  most 
likely  explanation  for  the  discrepancy  between  the  measured  and  predicted 
setup  profiles. 
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Figure  16.  Comparison  of  model-predicted  wave  decay  with 
and  without  bottom  friction  losses 


32.  As  noted  by  Bowen,  Inman,  and  Simmons  (1968),  the  measured  setdown 
is  nearly  uniform  in  the  region  between  incipient  breaking  and  the  point  where 
the  curl  of  the  plunging  breaker  touches  down.  It  is  interesting  that  even 
though  the  wave  height  is  decreasing  in  this  region,  the  momentum  flux  appar¬ 
ently  is  not.  Perhaps  this  is  because  (strictly  speaking)  no  energy  is 
dissipated  until  the  curl  touches  down  and  "white  water"  appears. 


SETDOWN/SETUP.TJfcm)  WAVE  HEIGHT,  H(cm) 


rr 


Figure  17.  Model-predicted  wave  decay  and  setdown/setup  as  compared  with 
laboratory  data  from  Bowen,  et  al.  (1968) 
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i.  Model-predicted  wave  decay  and  setdown/setup  as  compared  with 
laboratory  data  from  Bowen  (1969) 


PART  IV:  RESULTS  FOR  LARGE  SCALE  CONDITIONS 


33.  Applying  the  results  of  the  laboratory  calibration  to  near-prototype 
conditions  may  be  questionable  due  to  expected  scale  effects  in  the  turbulence 
associated  with  wave  breaking  as  previously  discussed.  However,  in  order  to 
demonstrate  use  of  the  model  for  waves  breaking  on  beach  profiles  containing 
bars  and  to  lend  some  validity  to  the  model  for  prototype  situations,  the 
model  was  run  for  large  scale  conditions.  A  "prototype  scale"  beach  profile 
generated  by  Saville  (1957)  in  the  Beach  Erosion  Board's  large  wave  tank  was 
utilized  and  is  displayed  in  Figure  19.  The  profile  is  characterized  by  three 
offshore  bar/trough  systems,  along  with  a  monotonic  section  in  the  nearshore 
region.  Test  conditions,  although  sketchy,  were  taken  from  the  laboratory 
notes  as 

Wave  period,  T  =  11.33  sec 

Wave  height  at  incipient  breaking,  =  1.78  m 

Location  of  outside  breaker  line  =  76.8  m  from  still-water 

line  on  beach 

Mean  sediment  diameter,  D  =  0.2  mm 

These  conditions  placed  the  breaker  line  outside  the  crest  of  the  first  bar 
at  a  still-water  depth  of  1.78  m.  The  setdown  in  mean  water  level  given  by 
Equation  36  is  -0.11  m,  and  the  ratio  of  wave  height  to  mean  water  depth  at 
incipient  breaking  is  1.07.  T  was  set  equal  to  0.4,  and,  because  the  seaward 
faces  of  the  bars  were  steep,  K  was  set  equal  to  0.2.  Following  the  proce¬ 
dure  described  by  Kamphuis  (1975)  and  assuming  the  bottom  was  not  rippled,  the 
bottom  friction  factor  f  was  found  to  be  roughly  0.005.  Initial  runs  showed 
that  bottom  friction  caused  decay  on  the  order  of  millimeters  only.  Appar¬ 
ently,  bottom  friction  plays  an  important  role  in  wave  decay  only  under  non¬ 
breaking  conditions  well  outside  the  surf  zone,  and  friction  was  again  left 
out  of  the  model  for  this  test.  Using  these  incipient  conditions,  the  trans¬ 
formation  of  wave  height  was  generated  in  a  stepwise  fashion  across  the  entire 
surf  zone  (Ax  =  1.5  m)  until  the  mean  water  depth  became  less  than  0.25  m. 

The  distribution  of  model-predicted  wave  height  and  setup  are  also  shown  in 
Figure  19.  Note  that  the  wave  reaches  the  stable  criterion  in  the  outermost 
trough  as  might  be  expected,  but  does  not  shoal  enough  on  the  second  bar  to 
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Figure  19.  Test  of  wave  transformation  model  at  prototype  scale  on 
large  wave  tank  profile  from  Saville  (1957) 


reach  Che  incipient  breaking  criterion.  Instead,  it  continues  past  the  second 
trough  (losing  magnitude  slightly  due  to  the  deeper  water)  until  the  third  bar 
is  encountered  where  breaking  starts  again  and  continues  until  the  shoreline 
is  reached. 

34.  The  results  of  application  of  the  model  under  large  scale  conditions 
seem  reasonably  valid,  at  least  in  a  qualitative  sense.  The  example  has 
demonstrated  the  ability  of  the  model  to  describe  wave  breaking  and  reforma¬ 
tion,  a  commonly  observed  process  on  natural  beaches.  The  predicted  wave 
decay  and  setup  profiles  are  continuous  and  well-behaved  until  the  mean  water 
depth  becomes  quite  small  (h  <  0.25  m) ,  where  a  swash  zone  model  might  be  more 
appropriate . 


PART  V:  CONCLUSIONS  AND  RECOMMENDATIONS 

35.  Based  on  laboratory  data  collected  by  Horikawa  and  Kuo  (1966),  the 
parameters  found  to  affect  most  the  decay  in  wave  height  due  to  breaking  in 
the  surf  zone  are  the  wave-height-to-water-depth  ratio  at  incipient  breaking 
and  the  beach  slope.  Wave-induced  setup  in  mean  water  level  plays  a  smaller 
but  non-trivial  role  in  governing  the  shape  of  the  wave  decay  profile,  espe¬ 
cially  near  the  still-water  line.  The  similarity  model  H  a  h'  commonly  used 
by  the  coastal  profession,  appears  to  be  reasonable  only  on  steep  beaches 
(1/20  to  1/30  at  laboratory  scale),  and  the  0.78  criterion  predicts  with  mar¬ 
ginal  accuracy  for  a  1/30  slope  only.  However,  the  model  developed  herein 
appears  to  qualitatively  and  quantitatively  describe  wave  transformation  in 
the  surf  zone  due  to  shoaling,  breaking,  and  reformation  over  a  wide  range  of 

beach  slopes  (1/80  to  1/12)  and  incipient  conditions  (0.63  ^  (H/h.)  ^  1.67). 

b 

Its  greatest  assets  are  its  simplicity  and  ease  of  application.  Although  it 
is  most  successful  on  profiles  of  monotonic  shape,  it  may  also  be  employed 
when  multiple  bar/trough  systems  are  present. 

36.  The  model  predicts  maximum  setup  values  with  reasonable  accuracy  for 
two  test  cases  presented  by  Bowen,  Inman,  and  Simmons  (1968)  and  Bowen  (1969); 
however,  it  does  not  describe  the  distribution  of  setup  across  the  surf  zone 
satisfactorily.  Apparently,  the  basic  assumption  that  onshore  radiation 
stress  can  be  described  by  using  linear  wave  theory  is  not  a  valid  one  under 
near-breaking  conditions.  More  work  is  required  in  this  area. 

37.  From  calculations  based  on  the  work  of  Kamphuis  (1975),  it  can  be 
concluded  that  bottom  friction  plays  a  negligible  role  in  wave  decay  in  the 
surf  zone  for  most  naturally  occurring  conditions  when  compared  with  the  ef¬ 
fects  of  breaking  and  shoaling.  Bottom  friction  could  be  significant  in  near¬ 
shore  regions  that  have  very  mild  slopes  or  rough  bottoms. 

38.  The  model  has  certain  limitations  which  might  restrict  its  use. 
Reflection  is  not  included,  and,  although  the  model  appears  valid  even  for  a 
1/12  slope,  one  must  be  careful  when  interpreting  results  for  runs  on  steep 
beaches.  An  elementary  scheme  utilizing  reflection  coefficients  might  be  in¬ 
cluded.  A  breaking  phenomenon  not  handled  well  by  the  model  is  wave  "trip¬ 
ping"  which  commonly  occurs  when  incipient  breaking  conditions  are  reached  at 
or  near  the  peak  of  a  bar  or  submerged  structure  and  the  wave  travels  onto  the 
trough  before  breaking  becomes  fully  developed.  The  model  might  be  improved 
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by  including  also  the  transfer  of  energy  to  higher  frequencies  as  studied  by 
Sawaragi  and  Iwata  (1974).  However,  the  cause  of  this  phenomenon  is  not  fully 
understood  at  present.  To  better  represent  natural  conditions,  another  logi¬ 
cal  improvement  would  be  an  application  to  random  waves,  perhaps  similar  to 
the  analysis  of  Battjes  and  Janssen  (1978)  or  Thornton  and  Guza  (1983). 
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APPENDIX  A:  ANALYTICAL  SOLUTIONS  FOR  THE  BREAKING  WAVE  MODEL 


Horizontal  Bottom 

1.  The  proposed  expression  (Equation  10)  for  predicting  wave  decay  due 
to  breaking  in  the  surf  zone  (wave-induced  setup  and  bottom  friction  not  in¬ 
cluded)  can  be  solved  analytically  for  the  flat  shelf  case  shown  in  Figure  lb. 
Starting  with  Equation  10* 


d(ECg) 

dx 


=  =$  [ECg  -  (ECg) g ] 


(Al) 


making  the  "stable  wave"  assumption 

H  =  IV  (A2) 

s 

and  applying  linear  shallow-water  wave  theory,  the  expression  takes  the  form 


1/2  dHf  _  zK  2,1/2 
1  dx  h'  (H  h 


-  r2h'5/2> 


(A3) 


Letting  F  =  H  and  rearranging  yield: 


dF  K_ 
dx  h' 


F  =  KT2h' 


(A4) 


This  first-order  ordinary  differential  equation  is  solved  by  summing  a  homo¬ 
geneous  solution  and  a  particular  solution,  or 


F 


F  +  F 
H  P 


The  homogeneous  solution  obtained  from 


+  K-  F 
dx  h'  H 


=  0 


(A5) 


(A6) 


*  For  convenience,  symbols  and  unusual  abbreviations  are  listed  and  defined 
in  the  Notation  (Appendix  B) . 


Al 


is  given  by 


-Kx/h' 


FH=  Ae 


where  A  is  an  arbitrary  coefficient.  Let  the  particular  solution  Fp  be  a 
constant  C  and  from  Equation  A4 

C  =  r2h'2  (A8) 


and  the  complete  solution  has  the  form 


F  =  Ae~Kx/h'+  r2h'2 


Applying  the  boundary  condition  at  the  breaker  line 


F  =  F  =  H.  at  x  =  0 
b  b 


the  coefficient  A  is  determined 


A  =  nr  - 


r2h'2 


and  the  solution  becomes 


H2  =  ^H2  -  r2h,2)e'Kx/h'+  r2h'‘ 


H  /  j  \  2  r2l  /-Kx\( 

P  =l([\hVb  "  rJ  exp  VhW) 


(A10) 


(All) 


(A12) 


(A13) 


Plane  Beach 


2.  For  the  case  of  a  plane  beach  described  by  the  expression 


h’ 


(A14) 


where  the  origin  is  again  at  the  breaker  line  and  x  directed  onshore. 
Equation  8  can  be  solved  analytically  again.  Following  the  same  procedure  as 
for  the  horizontal  bottom  case,  the  following  equation  is  developed 


£  +  ^  o  -  icr2h'3/2 

dx  h 


(A15) 


2  1/2 

where  G  =  H  h'  .  Converting  the  derivative  with  respect  to  distance  to  a 
derivative  with  respect  to  still-water  depth  by  utilizing  Equation  A14  yields 


dG 

dh' 


k _  g  =  -Rr2 

mh'  m 


.3/2 


(A16) 


Let  the  solution  to  the  homogeneous  equation 


dGH 

dhT 


G„  =  0 


(A17) 


have  the  form 


Ah' 


n 


(A18) 


where  A  and  n  are  arbitrary  constants.  Substituting  into  Equation  A17 
produces 


A  Inh 


,n-l 


h' 

m 


(A19 


In  order  for  the  homogeneous  solution  to  be  non-trivial,  n  =  K/m  and 


(A20) 


gh  =  Ah’ 


Let  the  particular  solution  have  the  form 


Gp  =  Bh’ 


(A21) 


where  B  and  n  are  arbitrary  constants.  From  Equation  A16 


B  (nh'n_1-  -  h,n_1)=  —  h,: 

\  m  /  m 


(A22) 


For  the  equality  to  hold  for  all  values  of  h'  ,  n  =  5/2  and  therefore 


(!-l) 


(A23) 


G  =  Ah 


jK/m  Kr 


»(!  -  I) 


(A24) 


The  coefficient  A  is  determined  by  applying  the  boundary  condition 


G  •  s  ■  "foI/2  at  h’  -  h; 


(A25) 


H^h^1/2  Kr2h£5/2 

A  =  ,  ,K/m  +  "75  K\,  ,K/m 

hb  m\2  -  m/b 


(A26) 


Substituting  into  Equation  A24  and  putting  the  result  into  dimensionless  form 


where 


[ (K/m)-l/2] 


(1  +  a) 


1/2 


(A27: 


(A28> 


Note  that  if  K/m  =  5/2  ,  the  solution  fails.  For  this  case  let  the  par¬ 
ticular  4 solution  have  the  form 


Gp  =  V (h* )h' 5/2 


(A29; 


and  from  Equation  A16  we  find 


V(h')  =  -  |  r2  In  h' 


(A30) 


and  the  total  solution  is  now 


G 


Ah 


,5/2 


5  r2,  ,5/2 
"  2  F  h 


)ln  h' 


Applying  the  same  boundary  condition  as  before  yields 


A 


in  h' 

b 


(A31) 


(A3  2) 


and  the  final  result  in  dimensionless  form  for  the  special  case  of  K/m  =  5/2 


APPENDIX  B:  NOTATION 


Cg  group  velocity 

E  wave  energy  per  unit  surface  area 

Eflux  enersy  flux 

E'  specific  energy  of  flow 

f  drag  coefficient  (bottom  friction  factor) 

FF  Froude  number 

g  acceleration  due  to  gravity,  constant 
h  mean  water  depth  including  setup 

h'  still-water  depth 

H  wave  height 

K  dimensionless  decay  coefficient,  wave  decay  factor 
m  beach  slope 

t  time 

T  wave  period 

u  horizontal  fluid  velocity 

x  horizontal  coordinate 

Y  unit  weight  of  water 

r  dimensionless  coefficient  (about  0.35  -  0.40),  stable  wave  factor 
Ax  size  of  numerical  step 

e  error  function 

n  setup 

p  mass  density  of  water 

Subscripts 
b 


s 


corresponds  to  incipient  breaking 
corresponds  to  stable  wave 


